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CHAPTER  I 


INTRODUCTION 

1.1  Use  of  the  Model 

DARCOM's  SESAME  (Selective  Stockage  for  Availability,  Multi-Echelon) 
model  is  of  the  multi-item,  multi-echelon  type.  It  determines  by  means  of 
mathematical  optimizing  techniques  how  many  of  each  component  to  stock  at 
each  type  stockage  point  in  the  supply  system,  taking  into  account  the 
potential  impact  of  each  backordered  component  on  system  down  time.  SESAME 
will  stock  to  achieve  any  given  weapon  system  target  availability  at  least 
cost.  Its  application  is  by  weapon  system,  primarily  to  "ERPSL"  weapon  sys¬ 
tems.  These  are  systems  for  which  stockage  based  on  routine  supply  rules, 
geared  to  overall  supply  performance,  is  inadequate  to  permit  achievement  of 
desired  weapon  system  operational  availabilities.  Typically,  these  systems 
have  low  densities,  have  reliable  components,  and  are  of  great  military 
importance. 

The  SESAME  model  is  used  in  both  budget  and  production  modes.  In 
budget  mode  its  primary  purpose  is  to  develop  curves  showing  the  relation¬ 
ship  between  target  operational  availabilities  and  necessary  inventory 
investment,  in  production  mode  it  produces  punch  cards  by  which  the 
stockage  quantities  it  computes  are  entered  into  the  Provisioning  Master 
Record,  a  Commodity  Ccmmand  Standard  System  (CCSS)  file.  CSSS  programs  use 
this  information  in  preparing  Essential  Repair  Parts  Stockage  Lists  (ERPSL 's) 
l  6r  the  field  units,  and  in  computing  buy  requirements  for  the  wholesale 
system.  The  ERPSL' s  give  stockage  quantities  by  unit. 

Part  support  is  only  one  factor  determining  operational  availability. 

The  SESAME  model  has  been  incorporated  into  a  life  cycle  cost  model  [16]  so 
that  all  the  relevant  factors  can  be  considered  at  one  time,  interactively. 
When  SESAME  is  used  in  budget  mode,  the  impact  on  operational  availability 
of  such  factors  as  system  reliability,  and  maintenance  policies  and  times, 
has  already  been  determined . 

1.2  Development  of  the  Model 

For  the  most  part  the  model  is  a  synthesis  of  already  existing  mathe¬ 
matical  techniques  incorporated  into  a  package  designed  to  maximize  user 
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convenience,  and  compatibility  with  CCSS.  This  package  was  designed  by 
the  DARCOM  Provisioning  Technical  Working  Group,  which  is  a  user  group 
chaired  by  DARCOM  HQ  and  incorporating  various  elements  of  organizational 
expertise;  e.g.,  IRO  has  been  the  primary  source  of  mathematical  expertise 
and  overall  technical  guidance.  The  model  was  actually  programmed  at  IRO; 

US  Army  Armament  Readiness  Command  made  a  number  of  program  enhancements, 
and  US  Army  Missile  Command  and  Army  Logistic^  Mi;  cement  Systems  Agency 
developed  programs  to  provide  Input  data. 

While  the  SESAME  model  is  currently  the  most  widely  used  model  for 
Army  ERPSL  systems,  it  does  not  address  various  complicating  factors  which 
may  be  significant  for  particular  weapon  systems.  This  is  by  design,  since 
it  was  felt  that  an  all  purpose  model  would  be  too  complex,  more  difficult 
to  use,  and  more  difficult  to  maintain  error  free.  Other  comparable 
models  are  available  [4,  9,  18],  drawing  on  much  of  the  same  mathematical 
content,  but  with  different  strengths  and  weaknesses  than  SESAME.  Guide¬ 
lines  are  being  developed  for  when  each  model  is  most  appropriate,  and  when 
it  may  be  necessary  to  tailor  SESAME  or  one  of  the  other  models  to  a  partic¬ 
ular  weapon  system. 

1 • 3  Purpose/Oryanlzation  of  Report 

This  report  is  intended  to  complement  an  extensive  user  guide  [6]  being 
written  by  a  sub-committee  of  the  DARCOM  Provisioning  Technical  Working 
Group.  This  report  discusses  the  mathematical  content  of  SESAME,  The 
chapters  on  assumptions  and  mathematical  overview  are  intended  for  the 
general  technical  reader.  The  other  chapters  are  intended  for  the  inventory 
modeller. 

There  is  still  a  third  related  report  which  provides  actual  computer 
program  documentation  for  the  SESAME  program  [7]. 
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CHAPTER  II 


ASSUMPTIONS 


Table  2.1  is  a  list  of  complicating  factors  which  are  relevant  to 
the  design  of  a  multi-echelon  model.  Checked  factors  are  incorporated 


into  the  current  version  of  SESAME.  Factors  with  a 


were  considered  in 


model  design  and  can  be  handled  in  an  approximate  fashion. 

SESAME  could  be  adapted  to  handle  just  about  any  of  the  other  factors, 
in  some  cases  exactly,  in  other  cases  approximately.  Since  each  new  factor 
complicates  the  model,  this  hasn't  been  considered  desirable;  further 
evolution  of  SESAME  is  expected.  Chapter  7  discusses  the  difficulties 
of  a  mathematical  nature,  if  any,  presented  by  each  of  the  untreated  factors. 

The  rest  of  this  chapter  is  a  gloss  on  the  less  obvious  complicating 
factors,  and  on  the  use  of  the  "<u . " 

Support  Structure.  By  units  differing,  we  refer  to  differences  among 
units  of  the  same  echelon,  e.g.,  in  order  and  ship  times.  A  non-vertical 
supply  structure  is  one  in  which  an  echelon  has  a  maintenance  but  not  a 
resupply  mission.  It  stocks  to  support  its  own  mission.  SESAME  assumes 
echelons  are  not  bypassed  so  times  are  additive:  If  ship  time  (given 
there  are  no  backorders)  from  Depot  to  Direct  Support  Unit  is  30  days,  and 
from  Direct  Support  Unit  to  Organization  is  5  days,  SESAME  assumes  that  if 
nothing  is  stocked  at  DSU,  then  ship  time  from  Depot  to  Organization  is  35 
days.  If  echelons  can  be  bypassed,  this  time  could  be  less. 

System  Structure.  An  application  is  a  particular  use  of  a  component  in 
a  Bystem.  The  pre-processor  for  SESAME  [17]  rolls  up  applications  so  that 
there  is  the  same  one  line  of  Information  about  each  component.  An  inden¬ 
ture  level  refers  to  the  hierarchical  role  of  a  component  In  a  system. 
Typically,  a  second  indenture  component /part  is  used  to  fix  a  first  inden¬ 
ture  component  which  is  used  to  fix  the  end  item.  SESAME  computes  stockage 
on  lower  indenture  parts  based  on  economic  considerations,  but  does  not 
explicitly  model  their  contribution  to  system  down  time.  Under  considera¬ 
tion  Is  a  proposal  to  compute  and  output  the  expected  average  delay  in 
getting  these  parts,  so  that  it  could  be  incorporated  into  the  repair 
times  of  the  reparable  components.  There  would  be  one  delay  estimate 
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made  for  each  echelon,  being  the  average  delay  for  all  components/parta 
used  at  that  echelon  to  repair  other  components. 

Failure  Process.  Failures  occur  at  random  times.  Some  component 

■># 

failures  may  induce  other  components  to  fail  so  that  failures  are  not 
independent  and  there  are  multiple  removals;  or,  because  of  mechanic's 
error  or  limitations  of  diagnostic  equipment,  it  may  be  likely  that  if 
component  A  fails,  component  B  will  also  be  removed  and  replaced.  SESAME 
allows  for  differences  between  the  number  of  system  failures  and  total 
number  of  component  failures/removals  in  its  approach  to  computing  opera¬ 
tional  availability.  All  component  failure  rates  used  in  SESAME  are  actually 
removal  rates. 

Maintenance  Policy.  The  same  component  may  not  always  be  removed,  or 
may  not  always  be  repaired  at  the  same  echelon.  Preventive  maintenance  is 
another  reason  there  are  removals  without  failures  and  is  therefore  treated 
by  SESAME  to  the  same  extent  that  non-independent  removals  are.  Under 
cannibalization,  an  inoperative  system  is  used  as  a  temporary  source  of 
components,  if  they  are  needed  and  otherwise  unavailable. 

Formally,  it  does  not  matter  if  the  time  to  repair  a  system  includes 
time  to  remove  the  failed  component — it  has  no  effect  on  any  calculations. 

However,  if  this  time  is  significant,  and  if  diagnosis  can  be  made  before 
removal,  an  adjustment  could  be  made  to  improve  results  (Chapter  7), 

On  many  systems,  test  equipment  narrows  source  of  a  problem  to  a  set  of 
components.  These  must  be  sequentially  replaced  until  the  problem  is  found. 

Resupply  Considerations.  If  an  organization  (ORG)  requests  a  component 
from  a  direct  support  unit  (DSU)  ,  and  the  DSU  is  temporarily  cut  of  stock, 
the  ORG  wait  depends  on  what  is  due-in  to  the  DSU.  Some  multi-echelon 
models  assume  the  ORG  request  will  be  delayed  the  full  order  and  ship  time 
from  the  depot  to  the  DSU.  The  issue  of  independence  of  successive  resupply 
times  is  treated  in  detail  in  Chapter  5.  Under  lateral  resupply,  a  unit's 
needs  may  be  filled  from  stock  held  by  other  units  at.  the  same  or  lower 
echelons. 

Under  real  time  distribution,  the  supplier  overrides  first  come,  first 
served,  depending  on  inventory  positions  of  customers  at  the  time  of  issue.  * 
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Usage /Environment.  In  general,  SESAME  is  a  steady  state  model,  which 
does  not  accommodate  changes  in  its  inputs  over  time;  it  will  compute  the 
budgetary  impact  of  phased  deployment. 

The  standard  definition  of  operational  availability  refers  to  per  cent 
of  calendar  time  the  system  is  operational. 

Relationship  to  Other  Problems.  Operational  readiness  float  is  a  pool 
of  systems  held  at  DSU  and  issued  when  a  system  is  returned  to  DSU  for  repair 
there. 
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TABLE  2.1:  COMPLICATING  FACTORS 


(1)  SUPPORT  STRUCTURE 
✓  (A)  MULTI-ECHELON 

/  (B)  EACH  ORG  (USING  UNIT)  SUPPORTS  MORE  THAN  ONE  OF  THE  WEAPON  SYSTEM 
t/  (C)  UNITS  DIFFER  BY  AREA 

(D)  UNITS  WITHIN  THE  SAME  AREA  DIFFER 

(E)  NON-VERTICAL  SUPPLY  STRUCTURE 

(F)  ECHELONS  CAN  BE  BYPASSED  TO  SAVE  TIME 

(2)  SYSTEM  STRUCTURE 

/  (A)  MULTIPLE  COMPONENTS 

'u  (B)  COMPONENT  WITH  MORE  THAN  ONE  APPLICATION 
%  (C)  MULTIPLE  INDENTURE  LEVELS 

(D)  REDUNDANCY 

(3)  FAILURE  T. TOC ESS 

/  (A)  FAILURES  ARE  STOCHASTIC  (RANDOM) 

^  (B)  NOT  ALL  FAILURES /REMOVALS  ARE  INDEPENDENT 

(C)  TIME  BETWEEN  FAILURES  IS  NOT  EXPONENTIAL 

(D)  MEAN  FAILURE  RATE  IS  NOT  KNOWN,  ONLY  ESTIMATED 

(4)  MAINTENANCE  POLICY 

/  (A)  COMPONENTS  ARE  NOT  ALWAYS  REMOVED  AT  SAME  ECHELON 
/  (B)  COMPONENTS  ARE  NOT  ALWAYS  REPAIRED  AT  SAME  ECHELON 
a,  (C)  THERE  IS  PREVENTIVE  MAINTENANCE 

(D)  CANNIBALIZATION  IS  UTILIZED 

(E)  DIAGNOSIS  CAN  OCCUR  BEFORE  REMOVAL 

(F)  FAILURE  IS  DIAGNOSED  BY  SEQUENTIAL  REPLACEMENT 

(5)  RESUPPLY  CONSIDERATIONS 

/  (A)  TIME  ON  BACKORDER  DEPENDS  ON  DUE-IN  POSITION  OF  THE  SUPPLIER 
/  (B)  SUCCESSIVE  RESUPPLY  TIMES  FOR  A  COMPONENT  ARE  NOT  INDEPENDENT 

(C)  USE  OF  EXPEDITED  RESUPPLY 

(D)  BATCH  ORDERING  OR  REPAIR 

(E)  LATERAL  RESUPPLY 

(F)  REAL  TIME  DISTRIBUTION 


TABLE  2.1  (CONT) 


(6)  USAGE/ENVIRONMENT 

y'  (A)  SYSTEM  OPERATIONAL  AVAILABILITY  MUST  BE  DETERMINED 
'v  (B)  DEPLOYMENT  IS  CHANGING 

(C)  USAGE  RATES  ARE  CHANGING 

(D)  NON-STANDARD  DEFINITIONS  OF  OPERATIONAL  AVAILABILITY 

(E)  STATES  OF  PARTIAL  DEGRADATION  MUST  BE  MODELLED 

(7)  RELATIONSHIP  TO  OTHER  PROBLEMS 

(A)  COMPONENTS  ARE  COMMON  TO  OTHER  SYSTEMS 

(B)  MOBILITY  OR  OTHER  CONSTRAINTS  0*»  COMPONENTS  STOCKED 

(C)  OPERATIONAL  READINESS  FLOAT  IS  UTILIZED 

(D)  REPAIR  CAPACITY  IS  LIMITED 


CHAPTER  III 


MATHEMATICAL  OVERVIEW 


3. 1  Basic  Approach 

The  objective  of  the  model  is: 


Minimize  Z  I  Stock  (I,J)  x  N(J)  x  Unit  Price  (I) 
I  J 


Subject  to  PNORS  <_  a 


where 


Stock  (I,J) 
N(J) 

Unit  Price  (I) 
PNORS 


a 


amount  of  item  I  stocked  at  an  echelon  J  unit 
number  of  stocking  units  at  echelon  J 
unit  price  of  item  I 

7.  of  time  system  is  down  due  to  unavailability 
of  a  component 
maximum  permiasable  PNORS 


The  PNORS  constraint  is  included  in  the  objective  function  by  restating 
the  problem  as: 

Minimize  Z  Z  Stock  (I,J)  x  N(J)  x  Unit  Price  (I) 

I  J 


+  Z  £  Expected  Backorders  (I,J)  x  RTD(I,J)  x  N(J)  x 
I  J  Backorder  Penalty  Cost (I) 


where 

Expected  Backorder  (I,J)  ■  expected  amount  of  item  I  backordered  at 

echelon  J 


RTD(I,J)  ■  replacement  task  distribution  percent. 

The  replacement  task  distribution  is  a  standard  Army  provisioning  element 
indicating  where  the  component  is  removed  and  replaced;  e.g.  RTD(I,1)  -  100% 
means  component  is  used  entirely  at  echelon  1  (ORG) . 

There  ara  now  three  questions: 

(a)  How  is  sytem  operational  availability  calculated  once  the 
model  is  solved?  Solving  the  model  means  finding  the  stockage  quantities 
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for  each  component  which  minimize  the  sum  of  backorder  and  investment  coats. 
Operational  availability  is  the  percent  of  calendar  time  the  system  is 
operational. 

(b)  How  are  backort’  *r  penalty  costs  determined  and  why  does  solving 
the  restated  problem  really  j'ive  a  solution  to  the  original  problem? 

(c)  How  is  the  optimum  solution  to  the  restated  problem  found? 

These  questions  are  discussed  non-mathematically  in  this  chapter,  and  then 
in  more  detail  in  subset  iert  chapters. 

3.2  System  Operational  Availability 

From  the  expected  backorders  for  a  component,  it  is  possible  mathema¬ 
tically  to  determine  the  "average  logistics  down  time"  at  user  level  for 
that  component,  i.e.,  the  average  time  to  get  that  component  when  it  is 
needed.  By  weighting  by  the  yearly  removal  rates  of  each  component,  it  is 
possible  to  determine  an  overall  average  logistics  down  time.  This  is 
combined  with  average  time  between  system  failure  and  with  down  time  while 
system  is  in  repair  to  compute  operational  availability. 

In  these  calculations  only  LRU's  are  considered.  LRU  was  originally 
an  Air  Force  term  meaning  line  replaceable  unit.  To  SESAME  users  it  is  an 
essential  component  whose  removal  and  replacement  restores  the  system  to  an 
operable  condition,  as  opposed  to  a  component  used  to  fix  a  higher  assembly 
which  itself  is  a  component.  A  carburetor  might  be  an  LRU  while  a  float 
valve  used  to  repair  the  carburetor  is  not  an  LRU. 

In  Chapter  II,  section  on  System  Structure,  treatment  of  non-LRU's  was 
discussed.  In  that  chapter  they  were  less  accurately  referred  to  aa  lower 
Indenture  components.  In  some  cases  a  second  indenture  component  could  be 
an  LRU  because  it  and  not  the  first  indenture  item  is  always  removed  and 
replaced  when  it  fails.  The  user  guide  provides  additional  information  on 
LRU  coding,  and  the  designation  can  be  made  by  an  automated  program. 

3.3  Backorder  Penalty  Costs 

Backorder  penalty  depends  on  whether  or  not  the  item  is  an  LRU.  The 

* 

backorder  penalty  for  all  LRU's  is  the  same  since  unavailability  of  any  LRU 
•k 

See  Chapter  VI  for  a  minor  qualification. 
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has  the  same  effect:  it  downs  a  system.  The  penalty  cost  is  called  the 
"CURPAR"  for  reasons  to  be  explained. 

When  the  restated  problem  is  solved  for  a  particular  CURPAR  there  Is, 
corresponding  to  the  solution,  a  total  inventory  investment  and  total  ex¬ 
pected  backorders.  Since  it  is  an  optimum  solution,  it  is  known  no  other 
pattern  of  stocking  costing  less  could  result  in  fewer  expected  backorders. 
We  have  a  least  cost  solution,  but  expected  backorders  may  not  correspond 
to  our  target  operational  availability.  If  the  backorders  are  too  high, 
we  raise  the  backorder  penalty  (CURPAR)  and  get  another  solution.  The  new 
solution  will  spend  more  on  stockage  since  each  backorder  avoided  now 


reduces  costs  by  a  greater  amount,  the  higher  CURPAR.  A  curve  is  developed 


Operational 

Availability 


Stockage  Cost 


,  X2,  X.j  correspond  to  solutions  found  with  successively  higher  CURPAR' s. 

From  the  curve,  the  CURFAR  needed  to  just  achieve  the  Operational 
Availability  is  found.  Since  SESAME  is  so  fast,  little  computer  time  is  lost 
in  generating  a  whole  curve.  The  CURPAR  selection  process  could  easily  have 
been  automated,  but  the  curve  itself  is  of  Interest  to  the  decision  maker. 

The  backorder  penalty  for  a  non-LRU  is  equated  to  the  cost  of  its  next 
higher  assembly.  The  rationale  for  this  is  that  if  the  average  number  of  back¬ 
orders  for  a  non-LRU  is  increased  by  one,  we  must  invest  in  at  least  one  addi¬ 
tional  next  higher  assembly  to  compensate  for  the  additional  assembly  lying 
unused  somewhere  awaiting  the  non-LRU  so  it  can  be  fixed. 

3 . A  Determining  Expected  Backorders 

This  is  highly  mathematical.  Backorders  depend  on  the  demand  rate  at 
user  level,  stockage  at  user  level,  repair  turn  around  time  if  the  item  is 
repairable  at  user  level,  and  order  and  ship  time  to  get  the  item  from  the 
user's  supplier  at  the  next  echelon.  The  order  and  ship  time  depends  not  only 
on  transportation  times,  but  on  whether  the  next  echelon  supplying  unit  is 
in  stock.  SESAME  actually  works  from  the  top  echelon  downward.  Suppose 
the  support  structure  consists  of  ORG,  DSU,  wholesale.  First  wholesale 
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performance  is  determined.  This  permits  us  to  calculate  the  order  nnd 
ship  times  the  DSU  will  experience.  Then  DSU  supply  performance  is  cal¬ 
culated  so  we  can  determine  the  order  and  ship  times  the  ORG  will  experience 
and  calculate  user  backorders. 

3.5  Finding  the  Optimvm  Solution 

The  objective  function  -  the  expression  for  what  is  being  minimized  - 
is  what  a  mathematician  calls  separable:  for  a  given  CURPAR,  the  solution 
for  each  component  can  be  found  separately  from  the  solutions  of  every  other 
component.  This  greatly  simplifies  things,  yet  tradeoff  between  components 
Is  still  accomplished.  We  illustrate  with  two  systems,  each  with  two 
components,  and  identical  input  data  except  for  the  unit  price  of  the  second 
component  (the  second  component  is  different  on  the  two  systems,  but  has  the 
same  failure  rates,  etc.). 

System  A  System  B 

Component  1  $150  $150 

Component  2  $100  $600 

Suppose  also  we  have  that  for  either  system: 

Alternative  Stockage  Quantities  to 
Achieve  Target  Availability 

Component  1  58 

Component  2  43 

i.e.  by  stocking  more  of  component  1  we  can  make  do  with  less  of  component  2. 

The  optimum  stockage  for  system  A  is  (5,4)  at  a  cost  of  5  x  $150  +  4 
x  $100  *  $1150.  For  system  B,  (5,4)  would  cost  $3150  while  (8,3)  costa 
$3000,  so  (8,3)  is  optimum. 

SESAME  would  reach  the  solution  as  follows: 


Stockage  System  A 

Stockage  System  B 

CURPAR  -  1 

3,2 

3,1 

CURPAR  -  2 

5,4 

5,2 

CURPAR  -  3 

- 

8,3 

5  of  component  1  and  4  of  component  2. 
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For  a  given  CURPAR,  SESAME  stocks  less  of  component  2  on  System  B  than  on 
System  A  because!  while  all  other  inputs  are  the  saute,  its  unit  price  is 
higher.  Therefore,  CURPAR-2  is  not  high  enough  to  pro luce  a  solution  for 
system  B  which  will  achieve  the  target  operational  availability.  Going  to 
CURPAR-3,  more  of  component  1  is  now  stocked  as  well  a*  more  of  component  2. 

For  any  given  component,  and  given  CURPAR,  the  least  cost  solution  is 
found  by  trying  all  possible  solutions  which  might  be  optimum  and  comparing 
costs.  This  is  feasible  because  the  computer  can  * ^iluate  each  solution 
very  quickly,  and  It  Is  possible  to  use  mathematics  co  rule  out  many 
solutions  as  not  being  optimum  without  actually  evaluating  them. 
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CHAPTER  IV 


CALCULATION  OF  OPERATIONAL  AVAILABILITY 


4.1  Notation 

OA:  "Operational  Availability";  hours  system  is  In  an  up  status  as 
a  per  cent  of  total  hours. 

SA:  "Supply  Availability";  per  cent  of  hours  syeietn  Is  not  down  due 
to  unavailability  of  a  component. 

IA:  "Inherent  Availability";  per  cent  of  hours  system  is  not  down  due 
to  need  to  remove  and  replace  an  (available)  component. 

EMF:  "Effective  Maintenance  F.n  w.r";  number  of  LRU  removals  per  end 
item  per  year. 

SFF:  "System  Failure  Frequency  number  of  times  system  goes  down  per 
year. 

LDT:  "Logistica  Down  Time";  average  time  to  get  an  LRU  when  it  is 

needed  to  repair  a  system. 

MTTR:  "Mean  Time  To  Repair";  time  to  repair  the  system  if  all  needed 
components  are  available. 

MCTBF:  "Mean  Calendar  Time  Between  Fa '.lures";  average  time  from  when 
system  goes  up  until  next  failure, 

"i":  Subscript  "i"  attached  to  any  of  the  above  denotes  component  i. 

4.2  Issues  Involved 

There  are  two  alternative  formulas  for  computing  operational  availability. 
To  rigorously  derive  either  formula  it  mutt  be  assumed  that  failures  are 
independent  and  there  are  no  multiple  removals  of  components  (see  Chapter  2) . 
What  we  label  the  "Case  A"  formula  also  assumes  cl.ut  while  the  system  is 
down  because  of  component  i,  no  other  comports  v  car.  fail.  The  "Case  B" 
formula  assumes  failures  of  other  component  ..a  eie  unaffected.  In  the  real 
world,  situations  exist  corresponding  to  hoi.h  Cane  A  and  Case  B  and  to 
mixtures . 

In  this  chapter  we  present  formulas  for  both  cases,  present  evidence 
that  choice  of  formula  does  not  have  much  effect  on  the  OA  computed,  and 
rationalize  choice  of  the  Cane  A  formula  as  likely  to  be  more  robust. 

15 


But, first  we  provide  some  background  information  on  the  relationship 
between  EKF  and  MCTBF. 

For  ease  of  exposition,  we  assume  all  system  failures  result  in 
component  removals-  If  they  do  not,  imagine  a  component  "*”  with  zero  LDT 
and  failure  rate  equal  to  the  rate  of  system  failures  without  removals. 

4.3  EMF  vs.  MCTBF:  A  Digression 

The  following  relationships  hold  for  both  Case  A  and  Case  E  for 
general  continuous  failure  distributions. 


(4.3.1) 


MCTBF  - 


MCTBF, 


(4.3.2) 


_0A_  .  0A  ,  _1_ 

MCTBF  L  MCTBF, 


Equation  (4.3.1)  states  that  time  between  failures  is  the  reciprocal 
of  failure  rate,  and  that  system  failure  rate  is  the  tium  of  component  failure 
rates,  i.e.  ,  the  (1/MCTBF^,) .  Equation  (4.3.2)  states  that  the  number  of 
system  failures  per  year  ie  the  fraction  of  a  year  a  system  is  up  divided 
by  the  mean  time  between  failures.  It  also  implies  (2nd  equality)  that  the 
proportion  of  system  failures  caused  by  component  i  is  the  relative  failure 
rate  of  component  i,  i.e.,  if  component  A's  failure  rate  is  twice  component 
B's,  it  will  cause  twice  as  many  system  failures  as  component  B.  This  is 
not  obvious. 

The  difference  between  Case  A  and  B  is  that  in  Case  A  all  removals 
correspond  to  system  failures  while  in  Case  B  component  failures  and  removals 
can  occur  when  the  system  is  already  down  (in  fact,  a  percent  equal  to 
100%-QA  occur  when  the  system  is  down) . 


(4.3.3a) 


EMF  -  SFF 


MCTBF 


Case  A 


(4.3.3b) 


SFF  -  (EMF) (OA) 


EMF  - 


OA  MCTBF 


Case  B 


Relationships  were  derived  by  Barlow  and  Proschan  [1]  and  Ross  [21]  as  cited 
in  Barlow  and  Hodes  [2]. 
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4.4  Calculation  of  OA,  SA:  Case  A 


Formulas  used  are: 


(4.4.1)  OA  -  (MCTBF) /(MCTBF  +  LDT  +  MTTR) 

(4.4.2)  SA  -  1  -  (LDT)/ (MCTBF  +  LDT  +  MTTR) 

MCTBF/ (MCTBF  +  LDT) 


Equation  (4.4.1)  is  well  known  (cf  Barlow  and  Proschan) .  Equation  (4.4.2) 
follows  from  the  same  arguments.  The  approximation  is  used  to  permit  calcu¬ 
lation  of  SA  independently  of  MTTR.  The  ratio  of  the  precise  value  to  the 
approximate  value  is,  after  algebra: 

(4.4.3)  Rat-o  -  1  +  (MCTBF) (MCTBF  +  LDT  +  MTTR) 

which  is  close  to  1  for  high  OA  systems;  e.g.,  if  MTTR/MCTBF  and  LDT/MCTBF 
are  each  5%,  ratio  is  about  1.0025. 

To  use  the  formulas,  the  LDT^  are  calculated  (Chapter  V).  Then 

(4.4.4)  LDT  -£(LDT1/MCTB?1)  i  l  (1/MCTBF  ) 

i.e.,  the  LDT^  are  weighted  by  failure  rates.  This  is  justified  by 
the  observations  made  concerning  the  meaning  of  the  second  equality  of  (4.3.2). 


4,5  Calculation  of  OA,  SA:  Case  B 


Formulas  are: 

(4.5.1)  OA  -  J 

(4.5.2)  SA  -  J 


MCTBF1  +  LDT1  +  MTTR^ 
MCTBFi 

MCTBFi  +  LDT1 


In  other  words,  each  component  is  up  a  given  percent  of  time,  Its 
status  (up  or  down)  is  independent  of  the  status  of  all  other  components, 
so  system  availability  is  the  product  of  the  component  availabilities.  The 
formulas  do  assume  a  component  does  not  fail  during  storage  or  installation. 
As  a  matter  of  theoretical  Interest,  a  correction  to  allow  for  failure  during 
installation  is  derived  in  Note  1. 
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In  the  exact  form,  use  of  equation  (4.5.1)  would  require  input  of  the 
HTTR^  to  SESAME.  An  approximation  eliminates  this  need: 


(4.5.3) 


OA  ~  (SA)  (IA) 


MCTBF, 

TT _ L 


MCTBF 

TT  _ i 


1  MCTBF1  +  LDTt  i  MCTBF^  +  MTTR^ 


Interestingly  the  ratio  of  the  precise  value  to  the  approximate  is 
exactly  the  same  as  found  in  (4.4.3). 

Equation  (4.5.2)  can  be  recast  in  a  very  different  looking  form  as 
suggested  by  Bernard  Price.  It  can  be  shown  (cf  Chapter  V)  that 


(4.5.4) 


TWBi/n  -  (LDTi)(EMFi) 


where  TWB^  is  expected  time  weighted  backorders  for  the  i  component  and 
n  is  the  number  of  systems  supported.  Viewing  item  i  as  a  system  unto 
itself,  so,  for  example,  SFF  ■  EMF^,  and  referencing  equation  (4.3.2), 


(4.5.5) 


EMF1  = 


(MCTBF1)/(MCTBFi  +  LDTi  +  MTTR) 
MCTB  - 


MCTBF t  +  LDT^  +  MTTI^ 


Hence,  combining  (4.5.4)  and  (4.5.5). 


(4.5.6) 


MCTBF^  +  LDT±  +  MTTR±  MCTBF^  +  LDT± 


and,  referencing  (4.5.2) 


(4.5.7)  SA  (1  -  TWB  /n) 

i 


This  form  of  the  Case  B  approach  is  extensively  used  by  Logistics 
Management  Institute  in  its  multi-echelon  work. 

4.6  Comparison  of  Formulas 

Some  sample  SESAME  runs  were  used  to  compare  the  Case  A  formula  and 
Case  B  formula.  These  runs  had  only  supply  information,  so  the  comparison 

I 

was  made  between  the  SA  formulas.  MCTBF  was  determined  by  (4.3.3b),  so 
what  the  comparison  shows  is  how  much  difference  use  of  (4.4,2)  versus 
(4.5.2)  makes  when  Case  B  assumptions  are  correct. 
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There  was  not  much  difference,  as  seen  in  the  table  below: 


System 

Generator 

EMF 

T,DT 

Case  B  Formula 

Case  A  Formula 

70.56 

.464  Days 

91.4% 

91.8% 

Missile 

96.39 

.141  Days 

96.3% 

96.4% 

Tank 

4.98 

8.29  Days 

89.3% 

89.8% 

Different 

Missile 

26.91 

.061  Days 

99.55% 

99.55% 

4 1 7  Choice  of  Formula 

The  advantage  of  the  Case  A  formula  is  that  it  permits  use  of  an  estimate 
of  MCTBF  which  is  not  derived  from  or  implicit  in  the  MCTBF^.  For  example, 
suppose  a  number  of  non-LRU  components  are  incorrectly  coded  as  LRU's.  This 
would  erroneously  increase  the  estimate  of  EMF,  but  the  estimate  of  MCTBF, 
which  might  originally  be  a  system  specification,  and  is  then  verified  by 
testing,  would  still  be  correct.  Using  the  Case  A  approach,  the  only  error 
from  a  miscoded  LRU  is  that  the  calculation  of  LDT  would  average  in  LDT^ 
not  corresponding  to  LRU's,  which  might  or  might  not  affect  the  average. 

Under  the  Case  B  approach,  availabilities  of  the  non-LRU '9  would  erroneously 
be  factored  in,  lowering  the  overall  estimate  of  OA. 

If  there  are  multiple  removals,  use  of  the  Case  A  approach  errs  in  that 
the  expected  time  to  get  n  components  is  larger  than  the  time  to  get  only 
one,  yet  Case  A  formula  would  input  an  LDT  based  on  the  latter.  Under 
some  circumstance  this  might  not  even  be  a  source  of  error.  Multiple 
removals  may  arise  when  components  are  replaced  on  a  trial  and  error  basis 
until  the  one  faulty  part  is  found.  If  an  unavailable  part  is  skipped  in 
the  search,  it  is  only  the  LDT  of  the  faulty  component  that  is  of  concern. 

In  any  event,  there  are  various  ways  to  estimate  the  LDT  to  get  n  components 
(cf  Chapter  7). 

The  Case  B  approach,  under  multiple  removals,  once  again  errs  in 
implicitly  assuming  that  each  removal  corresponds  to  a  system  failure 
with  probability  OA  (cf  A, 3. 3b). 
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CHAPTER  V 


CALCULATION  OF  EXPECTED  BACKORDERS 
5 . 1  Alternative  Approaches 

Expected  backorders  are  calculated  for  one  component  at  a  time.  There 
are  a  number  of  procedures  for  doing  this  including  METRIC  [22],  ACCLOGTROM 
[9],  TWOPT  and  Simple  Simon  [14].  Simple  Simon  is  most  precise,  while 
TWOPT  performs  best  among  the  other  three,  and  requires  significantly  less 
computer  time  and  storage  than  Simple  Simon.  It  is  likely  to  be  more 
adaptible  than  Simple  Simon  to  inclusion  of  additional  complicating  factors. 
Both  TWOPT  and  Simple  Simon  modules  are  available  to  users  of  SESAME. 

The  source  of  error  in  METRIC  is  its  assumption  that  the  order  and 
ship  times  experienced  are  Independent  from  requisition  to  requisition. 
Clearly,  if  you  learned  that  a  requisition  was  backordered  10  days,  this 
conveys  information  about  what  delay  a  requisition  placed,  say  one  day 
later,  will  experience.  The  independence  assumption  causes  METRIC  to  over¬ 
estimate  achievable  item  availability. 

The  source  of  error  in  ACCLOGTROM  i9  that  it  does  not  account  for  the 
role  of  due-ins  in  reducing  backorder  times.  This  was  discussed  in 
Chapter  II,  section  on  Resupply  Considerations, 

Simple  Simon  is  an  exact  model.  However,  it  assumes  Poisson  demand 
from  an  infinite  population,  i.e.,  it  does  not  account  for  a  reduction  in 
demand  when  a  system  goes  down.  Work  done  by  this  author,  such  as  that 
reported  in  [10],  indicates  this  is  a  minor  concern.  (See  Appendix  IV  of 
reference,  Table  III,  "Warm  Standby") 

Figures  1  thru  5  give  Notation  and  then  recursive  formulas  for  cal¬ 
culating  expected  backorders  by  each  of  the  approaches  mentioned.  Only 
Two  Point  will  be  discussed  here,  as  the  others  are  documented  elsewhere. 

The  ACCLOGTROM  approach,  as  shown  here,  is  actually  a  generalization  of 
the  approach  In  [9]. 

5 • 2  Explanation  of  TWOPT 

Notr.i-ron  is  in  Figure  5.1,  and  the  algebra  by  which  expected  backorderB 
nre  computed  is  In  Figure  5.2.  To  simplify  the  explanation,  it  Is  assumed 

20 


JKjfc 


4  * 


.  Zi^  W/T  ;s"r#  T^.-c  -: 


A?^V"^ai(6»?.liMiiin>a 


; 

i 

i 


i 

1' 


♦ 

I.  . 

I 

I 


all  removals  occur  at  user  level;  the  generalization  of  the  algebra  is 
straightforward.  We  will  also  assume  the  support  structure  consists  of  , 

ORG  and  DSU  echelons  and  a  depot  so  NECH  “3. 

Our  objective  is  to  calculate  TWB^,  To  do  so  we  use  equations  2a, 

2b  for  J  «*  1.  These  equations  require  in  addition  to  the  basic  data  - 
DENp  EMF,  P^,  TAT^,  PSUM^,  OST^  -  estimates  of  CDEL^  and  FILL^  where  CDEL^ 
is  the  average  delay  when  the  DSU  backorders  an  ORG  demand,  and  FILL^  is 
the  DSU  fill  rate.  To  get  these  estimates  we  must  use  equations  2a  thru 
2d  recursively  beginning  with  the  depot,  J  «=  3.  At  the  end  of  our  depot 
calculations  we  have  FILL^  and  DEI^,  from  equations  2c  and  2d  respectively. 

We  then  recycl  through  2a  thru  2d  for  J  -  2,  getting  FILL^  and  CDEL^  as 
required.  In  doing  our  depot  work,  we  need  FILL^  and  CDEL^,  but  these  are 
given  as  initial  conditions.  They  state  that  when  the  depot  orders  from  i 

the  manufacturer  there  is  a  fixed  lead  time,  with  no  chance  of  delay.  j 

Let  us  now  focus  on  2a  thru  2d.  The  due-in  to  a  stock  point,  from 
repair  and/or  from  the  next  echelon,  ia  represented  as  a  Poisson  process 
compounded  by  a  tw<.  point:  distribution  on  the  mean.  The  two  points  result 
from  a  simplified  representation  of  the  actual  continuous  distribution 
on  the  OST  experienced.  U^  corresponds  to  an  OST  equal  to  the  input  OST, 
which  assumes  no  delay,  and  corresponds  to  the  input  OST  augmented  by 
the  average  time  oil  backorder,  given  there  is  a  backorder. 

For  those  leaders  unfamiliar  with  Array  provisioning  terminology, 

] 

we  will  illustrate  the  computation  of  Uj^  for  a  DSU,  J  »  2.  Then  (DEt^) 

(EMF)  is  the  rate  of  removals  for  the  component  on  all  systems  supported  1 

by  the  set  of  ORG's  supported  by  a  DSU.  is  the  percent  of  these  removals 
which  result  in  unserviceables  being  returned  and  going  into  the  DSU  repair 
pipeline.  TAT  is  the  time  from  failure  until  the  unserviceable  has  teen  j 

repaired  and  returned  to  DSU  stock.  Hence  (DEN„) (EMF) (P,) (TAT?)  is  the  j 

i  t  L  j 

average  amount  due  in  from  repair.  j 

(1-PSUMj)  is  the  percent  of  all  removals  not  resulting  in  unservicables  j 

which  can  be  repaired  at  either  ORG  or  DSU,  i.e.,  the  percent  of  removals 
which  represent  draw-downs  on  DSU  assets  which  must  be  made  up  by  ordering 
from  the  depot.  (DEN,) (EMF) (1-PSUM2) (OST^)  is  therefore  expected  due-in 
to  the  DSU  from  depot. 
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Equation  (2b)  is  fiisc  justified  for  the  case  of  FILL^+^  ■  100%. 

Time  weighted  backorders  is  the  expected  number  of  backorders  at  a  random 
point  in  time.  The  probability  of  (X-S^)  such  backorders  is  the  probability 
that  total  due-in  is  X:  total  assets,  the  sum  of  on-hand  +  due-in  minus  back¬ 
orders  is  always  Sj,  so  if  due-in  is  X,  either: 

X  <  Sj.  .+  on  hand  is  -  X 
X  >  Sj  -+■  backorders  are  X  - 

Many  authors  have  shown  that  for  the  process  under  discussion  the  number 
due- in  from  repair  is  independent  of  the  number  due-in  from  procurement,  or 
in  this  case  the  next  echelon,  and  that  the  sum  of  due-in  is  Poisson  dis¬ 
tributed.  All  that  must  be  assumed  is  that  whether  a  given  removal  is  re¬ 
pairable  at  echelon  J  is  chosen  by  a  Bemouilli  trial  with  parameter  . 

Equation  (2b)  i3  only  approximately  correct  for  FILLJ+^,  j*  100%.  If 
we  are  at  a  random  point  in  time  t: 

CDEL _ OST _ 

a  b  t 

it  assumes  that  either  all  requisitions  initiated  in  the  interval  (a,b)  were 
backordered,  or  none  were  backordered.  In  this  3ense  it  overcorrects  for 
the  METRIC  assumptions  of  independence  between  requisitions.  However,  be¬ 
cause  it  uses  a  two-point  distribution  to  represent  a  continuous  OST  dis¬ 
tribution,  with  tails,  it  still  overestimates  availability. 

Derivation  of  (2c)  is  analogous  to  (2b),  reflecting  only  the  difference 
in  computing  fill  rate  and  time  weighted  backorders. 

Equation  (2d)  states  that  average  backorder  time  per  demand  is  (TWB)/ 
(Demands  Received)  so  that  delay,  given  there  is  a  delay,  is  (TWB/Demands 
Received) /FTlLj .  The  first  claim  is  the  well  known  "A  *  LW"  relationship 
most  simply  and  generally  proved  by  Stidham  [23].  The  second  claim  follows 
from: 

Ave  Delay  ■  (Delay  Given  There  is  a  Delay) (Probability  of  Delay) 

Note  that  since  LDT  is  average  delay,  we  must  multiply  by  FILL^  in  (3b) 
to  get  average  delay  back  from  CDEL^. 
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5.3  Evaluation  of  Alternative  Methods 

At  the  time  of  the  evaluation,  we  had  a  somewhat  limited  data  base 

consisting  of  301  Patriot  LRU's.  In  order  to  extend  the  generality  of  the 

analysis,  this  data  was  also  made  to  "resemble"  a  Howitzer  by  changing 

system  parameters.  Details  are  in  Table  5.1. 

To  derive  Table  5.2  each  of  the  approximate  methods  were  run  for  a 

range  of  CURPAR's.  For  each  solution  they  obtained,  the  stockage  determined 

* 

waa  also  evaluated  using  Simple  Simon.  Reported  in  Table  5.2  are  the  weapon 
system  supply  availabilities  projected  by  the  approximate  methods,  and  the 
exact  supply  availabilities  computed  by  Simple  Simon  for  the  same  stockage. 

In  the  Patriot  "MF  x  4"  run,  all  removals  rates  were  multiplied  by  4. 

Analysis  of  the  earlier  results  suggested  this  might  constitute  a  worst  case. 

Table  5.3  shows  actual  cost  comparisons  of  inventory  investment  needed 
to  achieve  various  levels  of  weapon  system  supply  availability  when  stockage 
is  determined  using  each  of  the  methods.  For  example,  on  the  Patriot,  it 
would  cost  $26.5  million  to  achieve  90%  weapon  system  supply  availability 
using  SESAME  with  the  METRIC  backorder  evaluation  modules  as  one  of  its 
subroutines.  Note  that  METRIC  would  project  higher  than  90%  for  $26.5 
million.  After  SESAME  with  METRIC  was  run,  Simple  Simon  was  used  to  calcu¬ 
late  the  solution  and  compute  true  availability. 

TWOPT  performed  very  well,  with  little  loss  in  accuracy  or  Increase 
in  cost  over  Simple  Simon.  Use  of  a  constraint  on  stockage  may  have  con¬ 
tributed  to  this.  In  Chapter  VI  two  constraints  on  optimization  used  in 
SESAME  are  explained.  In  this  research  the  SIP  lower  bound  constraint  was 

_ I  i.L_  _*.1 _ _ _ _ i. 

wocu  f  tuc  uuiicl  was  nut  • 
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i  W.  Karl  Kruse  designed  and  programmed  the  evaluator,  and  did  the  Simple 

’  Simon  runs  reported  In  Table  5.3  as  well  as  the  programming  for  those  runs. 
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All  variables  listed,  except  OSTj. 


and  DEN^ , 


vary  by  component. 


J  -  echelon;  J  *  1  is  user,  etc. 

Sj  -  stockage  (on  hand  +  on  order)  for  unit  at  echelon  J 

Uj  -  mean  pipeline  qty  (amount  due-in)  for  unit  at  echelon  J 

OSTj  -  order  and  ship  time,  to  echelon  J  when  echelon  J  +  1  is  in  stock 

t'ELj  -  average  increase  in  OST  because  of  backordering  at  echelon 
J  +  1 

CDEL  -  conditional  DEL  .  Delay  given  requisition  is  backordered. 

J  J 

TAT  -  repair  turn  around  time;  includes  time  to  evacuate  item,  if 
J 

necessary,  from  point  of  removal  to  echelon  J  as  well  as 
time  to  repair  it. 

Pj  -  Maint  Task  Distribution  fraction;  percent  of  total  removals 

which  are  repaired  at  echelon  J. 

J 

PSUM  -  j  P 

J  .1-1  J 

DENj  -  density  supported  by  each  unit  at  echelon  J 
EMF  -  effective  maintenance  factor  (demand  per  day  per  end  item) 

Dj  -  demand  on  unit  at  echelon  J 

TWBj  -  time  weighted  backorders  ■  expected  number  of  backorders  at 
a  random  point  in  time 

FlLLj  -  Fill  rate  provided  by  unit  of  echelon  J 

FILLj  -  1  -  FILLj 

p(x;u)  -  Poisson  density  for  parameter  u; 

NF-CH  -  number  of  stocking  echelons 
LDT  -  average  logistic  downtime  for  component 

FIGURE  5.1:  NOTATION 
2 1* 


Initial  Conditions 


FIGURE  5.3:  SIMPLE  SIMON 


FIGURE  5.4:  METRIC 


Initial  Conditions 


FIGURE  5.5;  ACCLOGTROM 


MF  IS  MAINTENANCE  FACTOR  OR  REMOVALS  PER  TOO  END  ITEMS  PER  YEAR. 
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CHAPTER  VI 


OPTIMIZATION 

6.1  Single  Item  Optimization 

The  procedure  discussed  here  was  developed  by  Kruse  [13]  for  a  two  echelon 
example  and  was  formally  generalized  to  more  than  two  echelons  by  Kotkin  [12]. 

The  objective  is  to  determine  upper  bounds  on  optimum  stockage  quantities, 
and  then  to  dynamically  lower  these  bounds  as  potentially  optimum  solutions 
are  evaluated. 

Notation 

Sj  «  stockage  at  echelon  J 

* 

TCj_^(S)  «  lowest  possible  sum  of  backorder  and  inventory  costs,  given 
Sj  ■  S,  where  inventory  cost  is  charged  only  for  stock  at 
echelons  1  thru  (J-l) . 

UP  ■  unit  price 

* 

Sj  -  optimum  stockage  at  echelon  J 
Procedure 

Procedure  is  illustrated  for  n^  echelon.  First  calculate  TC  ,(0) 

*  /  n—1 

and  TC  -.(<*>),  where  the  latter  Is  the  cost  when  an  echelon  n  stocker  is 
n-i 

always  In  stock,  and  the  former  is  the  cost  when  it  is  never  In  stock. 

Then 

(6.1.1)  (S*) (UP)  +  TC*  .(S*)  <  (0) (UP)  +  TC  .(0) 

n  n-i  n  —  n-l 

since  the  left  hand  eiee  is  the  total  cost  including  stockage  cost  at  echelon 

* 

n  if  S  “  S  ,  and  the  right  hand  side  is  the  total  cost  if  S  "0. 
n  n  n 

Hence,  using  (6.1.1) 

(6.1.2)  (3*) (UP)  +  TC  .(«)  <  TC  .(0) 

n  n-l  n-l 

A 

since  TC  (w)  <  TC  . (S  ) ;  i.e.,  as  upper  echelon  stock  is  raised,  delays 
n—l  n-l  n 
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which  by  (6.1.4)  provides  a  lower  bound  on  than  doeB  (6.1.3). 

Now  if  n  Is  the  upper  echelon,  we  try  all  values  for  S  until  we  reach 
*  n 
an  upper  bound  on  S  ,  keeping  track  of  the  lowest  cost  S  found  in  the  search. 

n  a  n 

For  each  value  of  we  find  TC  ^(S  )  by  solving  an  n-1  echelon  problem  in 

which  ship  times  are  a  function  of  Sn»  For  the  n-1  echelon  problem  we  use 

the  bounding  procedure  all  over  again,  so  that  a  problem  with  any  number  of 

echelons  can  be  solved  recursively. 

At  the  lowest  echelon  cost  is  a  convex  function  of  so  a  bounding 

procedure  is  not  required.  Also,  for  fixed  values  of  ST,  J  >  3,  Kruse  [13] 

A 

showed  is  a  (weakly)  monotonic  decreasing  function  of  S^.  This  fact  is 
used  to  save  computer  time. 

Implication  of  Users  Above  ORG.  A  percent  of  LRU  removals  may  occur  above 

ORG,  reflecting  weapon  system  evacuations  (refer  to  use  of  RTD(1,J)  in 

restated  problem  in  Chapter  III) .  The  bounding  procedure  must  still  work 

if  the  definition  of  TCt  .  (S)  is  changed  to  include  backorder  costs  incurred 
ttl 

at  n  echelon  if  S  is  set  to  S.  Backorder  costs  are  evaluated  as 
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(TWBj) (PCER) (CURPAR)  where 

TWBj  are  expected  backorders  at  nth  echelon 

PCER  is  percent  of  demands  on  echelon  accounted  for  by 
removals  there  rather  than  resupply  of  lower  echelons. 

6.2  Constraints  on  Optimization 

If  a  field  unit  experiences  a  frequency  of  demand  for  an  item  which 
equals  or  exceeds  its  stockage  criteria,  it  will  stock  that  item  in  an 
amount  at  least  equal  to  what  it  would  compute  using  routine  replenishment 
rules  regardless  ^f  what  SESAME  computes.  By  DoDI  4140.42,  this  entire 
amount  cannot  be  provisioned,  but  a  somewhat  lesser  amount  determined  by 
the  SIP  (Standard  Initial  Provisioning)  model  can  be.  For  items  which  will 
pass  stockage  criteria,  SIP  quantities  are  used  as  lower  bounds,  to  incorporate 
the  realities  of  field  stockage  into  the  SESAME  model.  For  most  ERPSL  systems 
this  can  have  little  impact  since  few  items  do  pass  stockage  criteria,  at 
least  at  lower  echelons,  and  SIP  stockage  is  limited.  SESAME  also  has  an 
option  to  use  a  lower  bound  which  approximates  post-provisioning  field  unit 
stockage  computations. 

Currently,  SESAME  is  not  used  to  compute  wholesale  quantities.  Instead 
retail  stockage  Is  based  on  initial  fill  rates  and  delays  actually  (.xperienced 
at  wholesale  level.  There  are  a  number  of  reasons  for  this  but  this  policy 
is  under  review. 

6. 3  Multi-Item  Optimization 

Case  A  Approach 

UBing  the  Case  A  approach  to  Weapon  System  availability,  our  ob¬ 
jective  is  to  minimize  inventory  investment  subject  to  a  bound  on  LDT.  We 
show  this  is  equivalent  to  constraining  on  the  sum  of  time  weighted  backorders, 
the  problems  we  actually  solve  (cf  Chapter  III) . 

Now,  recalling  equation  (4.4.4)  and  the  notation  of  that  Chapter 

(6.3.1)  LDT  -  E(LDTi/MCTBF1)/  E(1/MCTBF^) 
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Also,  by  (4.3.3a) 


(6.3.2) 


EMF  - 


OA 

MCTEF 


and  because  relative  removal  rates  are  proportional  to  relative  failure  rates, 
(cf  discussion  of  Equation  4.3.2)  , 

(EMF) (1/MCTBF ) 

(6‘3’3)  ”1  ~  --S-17MCTBF;  - 

Substituting  (6.3.3)  into  (6.3.1)  v 

(6.3.4)  LDT  -  E(LDTt)(EMFi) 

EMF 

Recalling  Chapter  IV, 

(6.3.5)  LDT^  -  TWB1/Demands  Received)  -  TWB^/(EMF^) (N) 

\ 

where  N  Is  number  of  weapon  systems  supported.  Hence,  by  (6.3.4)  and  (6.3.5) 

(6.3.6)  LDT  -  l  TWB1/(EM7) (N) 

Case  B  Approach 

Using  the  Case  B  Approach  to  operational  availability  we  would 
wiali  iu  minimize  (recall  Equation  4.5.7) 

TWB 

(6.3.7)  SA  -  it  (1  -  — 3 

i  w 

It  i 8  of  interest  how  a  solution  obtained  to  minimize  (6.3.7)  would  compare 
with  one  to  minimize  LDT. 

In  Appendix  2  we  show  that  the  solutions  must  be  close,  and  in  fact,  if 
in  solving  for  Item  i  CURPAR  were  adjusted  by  adding  UP^,  so  that  backorders 
were  costed  out  at  (CURPAR  +  UP^) ,  we  would  be  over  correcting  for  the 
multiplicative  objective  function. 
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This  adjustment  was  in  fact  made  in  SESAME,  at  a  time  when  the  Case  B 
approach  was  the  preferred  approach,  and  has  been  retained.  It  has  real 
impact  only  when  component  availability  targets  are  not  high,  and  helps 
moderate  extreme  differences  between  protection  afforded  high  cost  and  low 
cost  LRU's  in  that  situation.  This  seems  prudent  as  the  Case  B  approach 
might  be  more  correct  for  a  given  system,  and  if  there  is  any  cannibalization, 
extreme  differed ies  in  targets  is  not  desirable. 
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CHAPTER  VII 
EXTENSIONS 

7.1  Overview 

The  purpose  of  this  chapter  la  to  briefly  assess  the  mathematical 
difficulties  presented  by  those  complicating  factors  which  SESAME  does  not 
treat.  The  reader  is  referred  back  to  Chapter  II  and  its  discussion  of  what 
the  untreated  complicating  factors  are, 

As  in  any  such  chapter,  the  reader  proceeds  at  his  own  risk,  since  un¬ 
proved  and  untested  materiel  is  being  presented.  Furthermore,  the  whole  may 
be  greater  than  the  sum  of  the  parts;  difficulties  posed  by  treating  several 
factors  may  be  greater  than  implied  by  examining  each  factor  by  itself. 

7.2  Support  Structure 

Difficulties  presented  by  untreated  factors  are  not  of  a  theoretical 
nature. 

7.3  System  Structure 

Redundancy  of  Like  Items.  Underlying  "FSAME's  current  objective  function 
is  the  assumption  that  each  backorder  corresponds  to  another  down  system. 

Hence  if  pj  is  the  probability  of  j  backorders  (for  part  i) ,  and  each  system 
down  is  evaluated  at  a  coat  of  CURPAR, 

Cost  «  (CURPAR)  fl(j)  (p . ) ]  -  (CURPAR) [Ex(j)J 

i  2 

With  redundancy,  depending  on  the  redundancy  structure  and  the  number 
of  system  supported,  there  ie  some  non-linear  function  f(j)  which  maps  j 
into  the  expected  number  of  systems  down  given  j.  Calculation  of  itself 
is  affected  by  redundancy. 

Kaplan  [10]  discusses  the  computation  of  p ^ f Cj )  under  redundancy.  This 
expression  would  replace  Ex(j)  in  the  SESAME  model,  but  can  be  much  more 
difficult  to  compute  and  requires  knowledge  of  MTTR^  (time  to  remove  and 
replace  component  on  system) 

Redundancy  of  Unlike  Items.  In  this  case  the  system  can  operate,  for 
example,  if  either  component  i  or  component  j  is  up.  For  ease  of  exposition 
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assume  one  system  supported  per  ORG.  Using  the  Case  B  view  of  the  world 
(Chapter  IV),  which  ia  most  appropriate  in  this  context,  and  defining  and 
Aj  to  be  the  availabilities  of  components  i  and  j,  respectively,  system 
availability  is  A^  +  A^  -  A^A^  rather  than  A^A^ .  This  destroys  the  separ¬ 
ability  of  the  objective  function  (Chapter  III)  which  underlies  the  SESAME 
approach  to  multi-item  optimization. 

Presumably,  redundancy  of  unlike  items  would  Involve  (possibly  many) 
small  subsets  of  items,  i.e.  only  a  few  components  could  interchangeably 
do  any  one  function.  SESAME  could  approach  optimization  among  these  subsets 
by  iteration.  Continuing  our  example  of  two  components,  total  cost  for  the 
subset  would  be  evaluated  as  (Inventory  Investment)  +  (CURPAR) (A^  +  A^  -  A^A^) 
where  the  A, ,A  depend  on  the  p. ,p  ,  the  distributions  on  time  weighted 

^  J  ^  J  -I 

A-  **  *  1 _ J  1 _ _ -  1 _ 1  _ a _ J  «  4  A  .| _ J _ .  J _ 


backorders. 


Once  Aj  had  been  determined,  denote  this  A., ,  the  objective  for 


component  j  is  to  minimize  (Inventory  Investment  j)  +  (CURPAR)  (1-Ap  A^  . 

Once  the  A  resulting  from  this  minimization  is  found,  component  i  would  be 
J  2 

reevaluated  and  A^  determined.  More  generally  the  objective  is  to  minimize 

(Inventory  Investment)  +  (CURPAR)  (p, ,  p_...p  ).  The  iteration  procedure 

1  l  n 

would  stop  when  improvement  in  total  coat  dropped  below  a  threshhold.  No 
attempt  has  been  made  to  show  iteration  would  converge  to  the  optimum 
solutions. 

7.4  Failure  Process 

Time  Between  Failures  is  Not  Exponential.  The  simplest  generalization 
is  to  the  case  of  independent,  identically  distributed  interarrival  times, 
so  that  demand  constitutes  a  renewal  process.  If  N(t)  Is  the  number  of  re¬ 
newals  during  time  t,  the  distribution  for  M(t)  would  be  eubotituted  for  cue 
cumulative  Poisson  In  computing  backorders.  The  expression  for  N(t)  re¬ 
quired  is  that  for  the  equilibrium  renewal  process. 

Dependent  Removals .  Dependent  removals  might  result  in  demands  for 
more  than  one  unit  at  a  time.  This  has  been  discussed  in  the  literature  with 
justification  provided  for  use  of  Compound  Poisson  distributions  to  replace 
the  Poisson  in  the  evaluation  of  backorders,  within  the  context  of  the  METRIC 
(Chapter  V)  algorithm  [22].  Unfortunately,  we  require  not  the  average  time 
to  get  a  unit  of  stock,  but  the  average  time  to  get  all  units  demanded  in 
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one  requisition.  This  is  discussed  by  Kaplan  [11a]  for  single  echelon  work, 
and  is  readily  generalized  via  TWOPT. 

Chapter  IV  discusses  the  error  in  using  average  logistics  down  time 
when  more  than  one  component  type  is  required.  This  differs  from  requisition 
size  greater  than  one  since  the  times  to  get  each  cotoonent  required  are 
essentially  Independent.  Computation  of  the  time  to  get  n  components  is 
relatively  simple  for  reasonable  approximations  of  the  lead  time  process,  e.g. 
a  process  consisting  of  a  probability  of  zero  time,  and  a  uniform  or  ex¬ 
ponential  distribution  on  time  if  it  is  not  zero. 

Mean  Failure  Rate  is  Not  Known,  Only  Estimated.  Suppose  there  were 
two  possible  values  for  the  failure  rate  of  a  component,  with  known  prob¬ 
abilities.  This  can  easily  be  incorporated  into  the  optimization  by  evaluating 
all  stockage  levels  given  both  possible  demand  rates  and  weighting  backorders 
computed  by  the  probability  of  each  Tate.  This  is  not  equivalent  to  first  com¬ 
puting  one  expected  demand  rate  and  including  the  variability  due  to  unknown 
mean  in  total  variability  about  the  estimated  mean.  More  generally,  the 
Compound  Poisson  distribution  cannot  be  used  to  model  unknown  mean  at  each 
echelon  aa  is  sometimes  claimed  in  the  literature  for  this  would  not  capture 
the  correlation  between  mean  demand  at  each  echelon.  Since  SESAME  is  very 
fast,  an  entire  his* ogram  of  possible  true  means  can  be  accommodated. 

An  experiment  was  run  on  data  for  a  Missile  System  in  which  stockage 
levels  were  computed  assuming  that  the  failure  rate  was  known,  but  evaluated 

assuming  it  could  be  25%  higher  or  25%  lower  with  equal  probability.  This 

reduced  projected  system  supply  availability  from  89.5 %  to  86.4%  and  from 

99.4%  (with  a  higher  CURPAR)  to  98,6%.  Moreover,  three  high  failure  items 

had  to  be  excluded.  On  these  items  the  +  25%  experiment  reduced  availability 
precipitously;  for  the  lower.*  CURPAR  the  product  of  availabilities  for  the  three 
items  dropped  from  99.7%  to  91.4%. 

7.5  Maintenance  Policy 

Cannibalization.  Under  complete  cannibalization,  to  have  at  least  n 
systems  down  because  of  supply  unavailability  there  would  need  to  be  at 
least  n  backorders  for  one  of  the  components.  Without  cannibalization,  a 
total  of  n  backorders  on  n  different  components  could  produce  the  same  effect. 
In  SESAME  there  would  be  no  particular  difficulty  in  computing  the  entire 
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backorder  distril  tion  for  each  component .  Let 


F^(X)  '  probability  of  no  more  than  X  backorders  of  component  i 

?(X)  «  *  F  (X) 
i 

Then,  using  a  Case  B  type  approach  to  supply  availability,  and  a  well  known 
result  of  probability  theory  regarding  expectations, 

Prob  (No.  Systems  Down  >_ X)  »  1  -  F(X-l) 

m 

Ex(No.  Systems  Down)  «■  £  1  -  F(X-l) ,  m  ■  No.  of  Systems 

x*l 

Supply  Availability  ■  1  -  (Ex  No.  Systems  Down)/m 

The  approach  stated  does  not  consider  possibly  delays  due  to  removal 
of  the  component  being  cannibalized.  It  is  possible  to  model  this  delay 
accurately  if  it  is  assumed: 

a.  Component  is  only  cannibalized  when  needed  to  fix  a  system. 

b.  If  a  component  arrives  from  supply  it  is  installed  on  a  deadlined 
equipment,  even  if  it  is  not  sufficient  by  itself  tc  make  the  system  operable. 

c.  Once  cannibalization  is  begun,  system  down  time  equal  to  the  removal 
time  is  incurred. 

The  first  assumption  is  realistic,  the  second  may  be,  while  the  third 

is  not,  in  that  a  component  might  arrive  from  resupply  before  removal  was 

* 

completed.  Use  of  the  assumption  overstates  down  time. 

Given  the  assumptions,  the  proportion  of  the  year  a  system  is  down 
due  to  removal  of  component  1  may  be  calculated,  and  from  this  operational 
availability.  Let  OA,  and  IA  be  defined  as  in  Chapter  IV.  Also,  let  p^(x) 
be  the  probability  stock  on  hand  is  -  x,  and  let  R^  be  removal  time.  Delay 
due  to  removal  is  incurred  whenever  there  is  a  failure  of  component  i,  and 
stock  on  hand  for  i  Is  £  0  (negative  stock  denotes  backorders) ,  and  there  is 
a  system  from  which  the  component  can  be  cannibalized.  If  DR^  Is  yearly 
removal  time  per  end  item  and  there  are  EMFi  yearly  component  failures, 

DRi  -  (EKF  )(R  )  i  p  (x)[l  -  ^  F  (x) ] 

Jl»i  1 

*There  is  also  an  assumption  made,  of  little  consequence,  that  if  a  system 
fails  you  first  try  to  cannibalize  to  fix  it,  rather  than  cannibalize  it 
to  fix  another  system. 
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The  last  term  In  brackets  Is  the  probability  there  is  at  least  one  component 
with  at  least  X+l  backorders  so  that  there  are  (X+l)  systems  down  and 
component  1,  with  X  backorders,  can  be  cannibalized  if  a  demand  occurs. 

Then 

OA  (IA)  (1  -  Ex-J°^-Sy±tems  Down,  *  (  } 

n  i  i 

The  separability  property  of  the  objective  function  is  destroyed  by 
cannibalization,  making  optimization  difficult.  One  approach  would  be  to 
experiment  with  various  forms  of  non-linear  backorder  functions,  e.g.  set 
cost  -  (CURPARXTWBj)",  n  >  1. 

Diagnosis  Precedes  Removal.  Removal  time  must  be  considered  part  of 
"Mean  Time  to  Repair"  in  calculation  of  operational  availability.  Still,  if  a 
part  of  MTTR  precedes  logistics  dovn  time,  and  a  part  follows  it,  this  does 
not  impact  on  the  rationale  of  the  operational  availability  formulas.  How¬ 
ever,  if  there  is  built-in  diagnostic  equipment  so  diagnosis  precedes  removal, 
then  in  some  situations  logistics  down  time  and  time  to  repair  can  proceed 
concurrently. 

If  R^  is  removal  time  and  LDTi  is  logistics  down  time,  both  treated 
here  as  random  variables,  not  expectations,  then  adjusted  LDT^  is  max 
(0,  LDT^  -  R^) .  To  find  the  expectations  of  this  variable  requires  the 
distribution  on  adjusted  LDT^  (cf  discussion  on  Dependent  Removals) . 

7.6  Resupply  Considerations 

Expedited  Supply.  Suppose  that  d  days  between  DSU  and  ORG  can  be  saved 
by  expedited  resupply.  This  information  might  be  used  in  a  number  of  ways 
to  adjust  logistics  down  time  downward. 

a.  If  there  la  no  ORG  stockage,  subtract  d  days. 

b.  As  adjustment  a,  but  also  subtract  (d) (1-0RG-FILL)  (1-DSU-FILL)  if 
ORG  does  stock.  "FILL"  refers  to  fill  rate.  The  rationale  is  that  if  the 
ORG  is  out  of  stock,  and  there  is  not  due-in  from  the  DSU  already  enroute, 
i.e.  the  DSU  is  out  of  stock,  then  d  days  will  be  saved. 

c.  As  adjustment  b,  but  also  subtract  (1-0RG-FILL) (DSU-FILL) (f (d ,t) ) 
where  t  is  routine  resupply  time  between  DSU  and  ORG  and  f  is  some  function. 
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The  assumption  is  that  even  if  stock  is  enroute  to  the  ORG,  if  it  wi)l  not 
routinely  arrive  within  d  days,  expediting  can  help.  To  determine  f,  it 
might  be  assumed  time  to  receipt  of  the  due-in  is  uniformly  distributed 
between  0  and  OST,  or,  considering  the  possibility  of  multiple  due- ins,  a 
Beta  distribution  might  be  used  (if  there  are  n  Poisson  arrivals  in  an 

4.1. 

interval,  time  to  i  is  Beta  distributed). 

It  is  possible  to  incorporate  expedited  supply  directly  into  the 
optimization  routine. 

Batch  Ordering.  A  straight  forward  approach  would  be  to  compute  the 
economic  order  quantity  first,  and  then  optimize  on  the  total  stockage 
objective,  i.e.,  reorder  point  plus  order  quantity.  At  upper  echelons  the 
demand  process  would  be  Compound  Poisson,  to  account  for  demands  greater 
than  1  unit.  Problems  are: 

a.  Impact  on  validity  of  Two-Pt  approximation  approach.  Simple  Simon 
could  handle  exactly  a  two-echelon  problem  with  only  upper  echelon  using 
EOQ. 

b.  Demands  on  upper  echelon  are  correlated  if  lower  echelon  follow 
batch  ordering  policies. 

Muckstadt  [19]  has  developed  an  exact  solution  which  is  computionally 
tedious  and  has  never  been  implemented.  For  larger  number  of  customers  he 
suggests  ignoring  correlations.  Dumraeyer  and  Schwarz  [8]  propose  a  heu¬ 
ristic  which  they  tested  by  simulation. 

Lateral  Resupply.  Problems  could  be  horrendous. 

7. 7  Usage /Environment 

Barziiy  and  Gross  [ 3 ]  have  done  a  study  on  coping  with  changing  environ¬ 
ments  in  a  single  echelon  problem  with  limited  repair  capacity.  Bein  [4] 
has  considered  degradation.  While  different  definitions  of  operational  avail¬ 
ability  must  be  treated  on  a  case  by  case  basis,  there  is  an  entire  class  of 
weapon  systems  which  are  sortie  oriented,  so  that  the  availability  measure  of 
interest  is  availability  during  sortie  time.  At  one  extreme  such  problems  may 
degenerate  to  single  echelon  "fly-away"  kit  problems  in  which  resupply  is  not 
much  of  a  consideration.  This  occurs  if  time  between  sorties  is  large  so 
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that  resupply  can  always  occur,  while  sortie  time  Itself  is  small.  At  the 
other  extreme  this  procedure  might  he  used: 

a.  Compute  logistics  down  time  as  currently. 

b.  Based  on  an  approximate  distribution  on  down  time,  on  sortie  time, 
and  on  time  between  sorties,  compute  the  percent  of  down  time  which  occuro 
during  and  not  between  sorties. 

This  approach  Ignores  the  impact  of  sorties  on  the  failure  process. 

7.8  Relationship  to  Other  Problems 

Common  Items.  Suppose  there  were  two  systems  with  LRU's  distributed: 


One  approach  would  optimize  all  three  set3  individually,  over  ranges 
of  curve  parameters,  and  then  search  for  optimum  combinations  of  curve  para¬ 
meters.  Note  that  in  computing  expected  logistic  down  time,  the  contribution 
of  "common"  would  vary  between  weapon  system  A  and  B,  depending  on  the  non- 
common  LRU’s  as  well  as  the  failure  frequency  on  A  and  on  B  of  each  common 
LRU. 

A a  the  number  of  systems  grows,  the  number  of  sets  to  be  defined  grows 
geometrically,  e.g.  on  three  systems  there  are  parts  common  to  A  and  B,  B 
and  C,  A  and  B  and  C.  In  principle  it  should  be  feasible  to  cope  with  a 
great  number  of  sets,  but  a  simpler  approach  might  be  desirable.  Such  an 
approach  is  to  first  optimize  for  A,  but  also  compute  the  contribution  to 
down  time  on  other  systems  of  the  LRU’s  on  A  common  to  them.  Next  B  is 
optimized,  and  so  on.  In  doing  the  optimization  on  A,  all  backorders  are 
costed  out  at  the  CURPAR  being  used  for  A,  even  though  some  of  the  backorders 
impact  on  other  systems. 

Mobility  or  Other  Constraints.  Suppose  there  were  a  limit  on  the  number 
of  components  to  be  stocked  at  a  unit.  Using  the  Lagrangian  approach,  a  cost 
would  be  included  in  the  optimization  for  each  component  stocked,  and  this 
cost  would  be  raised  to  achieve  desired  stockage  list  size.  Since  now  there 
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are  two  Lagranglans,  pertaining  to  list  size  and  backorder  cost  (CURPAR) , 
the  search  process  for  correct  values  would  probably  need  to  be  automated. 

Operational  Readiness  Float.  This  is  treated  by  Kaplan  [11].  It  is  not 
particularly  difficult. 


jacts  on  Simple  Simon /TWOPT 


The  discussion  of  batch  ordering  raised  the  issue  of  how  the  factor 
impacted  on  Simple  Simon  or  the  TWOPT  approximation.  Simple  Simon  assumes 
Poisson  demand,  so  a  number  of  the  other  factors  also  bears  on  this.  Robust¬ 
ness  of  TWOPT  has  never  been  tested. 
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NOTE  1 


1  f _ 

p  f+r 


f  (x)  (f+r)e~(f'''r)dx 
o 


1  f_  1 _  1 

p  f+r  f+r  “  f+r 


When  repair  times  are  expo-r.ential ,  the  fact  that  repair  was  completed 
before  a  failure  occurred,  tells  us  something  about  hov  long  the  repair 
took.  If 

MTTRC  h  Ex  (Repair  Time [Repair  occurs  before  failure) 

Then,  by  a  derivation  analogous  to  derivation  of  t, 


MTTRC  - 


CD 

,  -fx  -rx, 

/  x  r  e  e  dx 

o _ 

1-p 


1 

f+r 


This  equals  t.  In  other  words,  something  will  happen  in  average  time  of 
(f+r) .  With  probability  p  it  is  a  failure  and  probability  (1-p)  it  is  a 
repair. 

TMTTR  -  (1-p) (MTTRC)  +  p[t  +  LDT  +  TMTTR] 
and  by  algebra 

TMTTR  -  MTTRC  +  (r2-)  [t  +  LDT] 

1"P 

■  ffe  »  +  f>  +  7  «-D« 

’  t  +  I  <LDT) 
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NOTE  2 


SOLVING  MULTIPLICATIVE  OBJECTIVE  FUNCTION 


Notation: 

TC  :  Total  Cost,  Inventory  plus  backorder 
Si  :  Total  stock  (for  item  1). 

UP^  :  Unit  price 
X  :  Lagranglan 

N  :  Number  of  weapon  systems  supported 

:  Expected  time  weighted  backorders  at  user  level 
SA  :  Supply  Availability 
Discussion: 

Our  objective  is: 


(1)  Minimize:  TC  -  r  S  UP  +  (XN)[1  -  *(1  -  -i)] 

i  i  " 

where  the  second  term  on  the  right  hand  aide  expresses  expected  number  of 
systems  down  given  the  Case  B  approach  to  operational  availability  (recall 
Equation  4.5.7) . 

Differentiating  with  respect  to 


(2)  9rc  <?  Bi 

as.  +  AV?s, 


where 


i-B  B 

(3)  At  -  it  — i  -  SA/(1  -  -~) 
i  ^  N  N 


Now  suppose  we  had  the  optimum  solution,  S. 


objective  as 


Then  if  we  restated  the 
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(4) 


Minimize:  TC  -  l  SiUP1  +  XAj^ 

* 

and  treated  Aj-  as  a  constant,  found  from  (3)  using  the  SA  determined  by  , 
we  would  find  that  the  same  minimized  (4)  as  had  minimized  (1) ,  This  is 
clear  since  the  derivative  of  (4)  with  respect  to  is  (2) . 

We  can  rewrite  (4)  as 

(5)  Minimize  TC  -  l  +  X^ 

where 


(6)  \±  -  XAj- 


Now,  we  will  show  later  that  given  an  optimum  solution  of  (5) , 


(7) 


Realistically,  UP^/X.^  will  not  exceed  1.  Using  this  fast  and  (6),  (3)  and  (7); 


(8)  (X)  (SA)  <  X± 


or,  by  algebra 


(X) (SA)  (X) (SA) 

1  -  Bi  (1  -  UP1/X1) 

~N 


/o\ 


/ 


\  < 

>'»/  > — * 


fX) (SA)  +  UP. 

1 


To  summarize,  we  have  found  that  solving  (1)  is  equivalant  to  solving 


(10)  Minirlze  TC  -  I  StUPt  +  x^ 
X'  <  \±  <  X'  +  upt 
X'  -  (X)  (SA) 
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The  significance  of  (10)  is  that  if  we  minimized  (10)  for  a  range  of 
X',  evaluated  Supply  Availability  for  each  solution  by  the  Case  B  approach, 
and  drew  the  availability  versus  investment  curve,  our  curve  and  accompanying 
stockage  lists  would  be  the  same  as  minimizing  (1) . 

Comparing  (10)  to  SESAME's  objective  function,  the  only  difference  is 
that  X  should  be  increased  by  a  maximum  of  UP^  when  item  i  is  being  solved. 
This  will  normally  have  very  little  impact  as  UP^  <  <  X  . 

Proof  that  B. /N  <  UP/X^ 

Let  be  expected  stock  unavailability,  defined  here  as  the  probability 
there  is  at  least  one  backorder  at  a  random  point  in  time. 

Then 


(11) 


Ui^ 


B 


i 

N 


For  example,  imagine  that  whenever  there  was  at  least  one  backorder,  there 
were  N.  This  is  a  worst  case.  In  this  case  «  B^/N  (e.g.  if  U ^  *  102, 

B1  -  (.10) (N)). 

Now  in  the  literature  there  is  a  well  known  result  which  characterizes 
when  it  corresponds  to  the  optimum  stockage  solution  in  a  single  echelon 
problem  with  Poisson  demand: 

UP 

<  Cost  Per  Unit  Asset  Per  Unit  Time  m  _ i 

^  i  Cost  Per  Unit  Backorder  Per  Unit  Time  ^ 

The  same  relationship  is  directly  extended  to  multi-echelons,  because  once 
the  optimum  solution  for  echelons  2  thru  NECH  is  fixed,  we  are  left  with  a 
classical  one  echelon  problem,  at  least  when  using  METRIC  or  TWOPT 
formulation. 

Combining  (11)  and  (12)  provides  the  desired  results. 
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